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Abstract 

We analyse the thermal evolution of the Universe in the Inert Doublet Model for three 

^, regions of Dark Matter masses: low mass (4 — 8) GeV, medium mass (30 — 80) GeV and high 

(^ mass (500 — 1000) GeV. Those three regions of DM mass exhibit different behaviour, both in 

CN the possible types of evolution and in the energy relic density values. In this analysis we use 

'^ the masses of the scalar particles as the input parameters to constrain the two self-couplings 

t— 5 between neutral scalars: A345,A2. These couplings are used to construct the parameter space 

T-H where different types of evolution may be presented. We discuss the influence of the scalar 

^~~^ masses on the type of the evolution. We also discuss the importance of the A2 self-coupling. 

I— I We argue, that the astrophysical data along with the positivity constraints simultaneously 

i-^ constrain both A345 and A2 self-couplings. 

^ 1 Introduction 

,__! Inert Doublet Model (IDM) is one of the widely discussed extensions of the Standard Model 
> (SM) that may provide the Dark Matter (DM) candidate P |2]. It is a Z2 symmetric 2HDM 

Qv with a suitable set of parameters. In this model two scalar doublets are introduced. $5 is the 

0\ Higgs doublet responsible for the electroweak symmetry breaking and masses of fermions and 
gauge bosons as in the SM. This doublet provides the longitudinal degrees of freedom of the 
gauge bosons and the SM-like Higgs boson particle hs- The other doublet, ^d, does not receive 
vacuum expectation value (v.e.v.) and does not couple to fermions. All components of this 

^ doublet are realized as the massive scalar D-particles: two charged D^ and two neutral D^ and 
K*- Da- The stability of the lightest of these particles, which is a candidate for the DM particle, 

S^ originates from the discrete symmetry of Z2 type, called below Z)-symmetry: 

D: $5 -^ $5, $B A -$/), SM fields A SM fields. (1) 

We assume that the current vacuum state of the Universe is described by IDM with the DM 
candidate Dh, so that: 

Md±,Md^>Md„. (2) 

In this paper we discuss the evolution of the Universe after infiation, following the approach 
presented in P, IH El El E] ■ In this approach the scalar potential changes due to the temperature 
corrections ~ T^ and in the past, for T 7^ 0, the vacuum structure could have different properties 
than the present vacuum state at T = 0. Sequences of the different vacua (represented on the 
phase diagrams by rays) describe different types of evolution of the Universe. 
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We consider the possible rays that can be reahzed in three regions of the DM mass: low 
mass (4 — 8) GeV, medium mass (30 — 80) GeV and high mass (500 — 1000) GeV. We use the 
estimation of the energy relic density of DM to show that for different mass regions different 
types of evolution are consistent with the WMAP constraints for the energy relic density of DM, 

In [6, 7j we focused on (/ii, ^2) phase diagram defined by the quadratic and two quartic parani- 
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eters of the scalar potential V12 
in a following way: 
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This phase diagram allowed us to present the possible types of evolution in easy way in form 
of rays. Each of the considered evolutions can be realized only if certain conditions for /ii, fi2 are 
satisfied. In the numerical analysis in [7] we have fixed the scalar masses, while the values of the 
scalar self-couplings A2 and A345 were chosen to fulfil the constraints for each ray. In this work 
we focus on those couplings, mainly on the A2 parameter, which is a quartic coupling for the DM 
particles, usually neglected in analysis of QnMh"^- However, this parameter has an important 
impact on the evolution and the history of the Universe. 

The content of this paper is as follows. In section 2 we present the main properties of 
IDM. In section 3 we discuss in details the DM relic density constrains in three regions of DM 
mass. In section 4 we introduce the (A345,A2) plane, which is useful in the discussion of the 
parameter regions for different rays, at the fixed scalar masses. Here we use the benchmark 
points to illustrate three regions of DM mass and present the connection between the type of 
evolution and the values of fi/jjv/^^- Section 5 contains a discussion on influence of different mass 
parameters on the (A345, A2) regions along with the relic density values. 



2 Model properties 
2.1 IDM 

Lagrangian We consider an electroweak symmetry breaking (EWSB) via the Brout-Englert- 
Higgs-Kibble (BEHK) mechanism described by the Lagrangian 



C = C'J' + CH + CY{i,f,<^, 



c 



H 



T~V. 



(4) 



Here, C^¥ describes the SM interaction of gauge bosons and fermions. 

The Higgs scalar Lagrangian Ch describes the interaction of two scalar doublets $5 and 
$£1 with the standard kinetic term T and the scalar potential V. The potential V, which can 
describe IDM, is invariant under D-symmetry: 



y = -: 



m?i$^$5 + "^22*l?*D 



+ 



1 



+A3 (<^l<^s) (^^d'^d) + A4 ($i$D) (^l,"^^ 



Ai (^$t;$s) + A2 (^$l3$D 



+ 1'' 



(5) 



All parameters of (tsl) are real and one can fix A5 < without a loss of generality. 

In order to have a stable vacuum we impose positivity constraints in the following form: 

Ai>0, A2>0, i?+l>0, 
A345 = A3 + A4 + A5, R = A345/V A1A2. 



(6) 

(7) 



They assure that the potential is bounded from below and the extreniuni with the lowest energy 
will be the vacuum (the global minimum of the potential). 

D-symmetric potential has 7 independent real parameters: ml^, ^^22' -^i-s- After EWSB 
those parameters can be expressed by the non-zero v.e.v, four scalar masses and two self-couplings 
between scalars. 

Cy describes the Yukawa interaction of fermions ipf with only one scalar doublet $5. It has 
the same form as in the SM with the change $ — )■ $5 (Model I for Yukawa interaction). Cy 
respects D-symmetry in any order of the perturbation theory. 

Inert vacuum We assume that the vacuum state of the potential (tsl) is given by the inert 
state (denoted by /i). In this case only $5 acquires the non-zero vacuum expectation value, 
which as it follows from extremum condition is equal to: 

v'' = mlJXi. (8) 

For Ji to be the vacuum following conditions should be satisfied |6]: 

/ii > for any i?, yUi > yU2 for i? > 1, i?yUi > ^2 for \R\ < 1. (9) 

The inert vacuum state is invariant under the /^-transformation just as the whole basic 
Lagrangian (|4]). There exist four dark scalar particles Dh^ Da,D'^, which are D-odd, and the 
Higgs particle hs, which interacts with the fermions and gauge bosons just as the Higgs boson 
in the SM. hs and the SM fields are D-even. In the inert vacuum the D-parity is conserved, and 
due to this fact the lightest D-odd particle is stable, being a good DM candidate. 

The masses of the physical fields hs, Dh, Da and D^ can be used to express the parameters 
of V after EWSB. Those relations are given by: 

Ml = W=mj„ Ml, = ^^"' : "^^^ , 

A A A +A ^^°^ 

Ml^=Ml, + ^v\ Ml^=Ml, + ^v^. 

Other crucial parameters after EWSB are the couplings between the scalars. To complete 
the set of parameters we use two of them, A345 and A2. The A345 is a coupling between SM-like 
Higgs hs and DM candidate Dh'- D^Duhs and DuDiihshs- The A2 coupling is proportional 
to a quartic self-coupling among D-particles, e.g. D^DhDhDh- The self-coupling that governs 
the charged scalars' interactions: D^D^hs and D^D~hshs is proportional to A3. 

2.2 Collider constraints 

Various theoretical and experimental constraints apply to the IDM (see e.g. [SI [HI [IHl IHl [121 [131 
[HI [ISl [ISl [HI [IB])- Below we list important existing constraints on couplings and masses both 
for the Higgs particle and the dark scalars. 

Constraints on self-couplings The positivity constraints are imposed directly on quartic 
parameters in the potential. If we want to assure the perturbativity of the theory, the self- 
couplings A's (i.e. A2,A3,A345) cannot be large. The bound (called perturbativity constraint) is 
set typically to 

|A|<47r. (11) 



Electroweak precision tests EWPT constrain strongly physics beyond SM. In IDM tliey 
limit the allowed values of masses both for the Higgs particle hs and the dark scalars. For IDM 
both light and heavy Higgs particle is allowed by EWPT j2]. Constraints for the dark scalars 
can be conveniently expressed by the limits for the mass splittings: 

6a = Md^ - Md„ , 6± = Md± - Md„ ■ (12) 

In [11| it was obtained that for a light Higgs boson, the allowed region corresponds to 6± ~ 6a_ 
with mass splittings that could be large. For heavy SM Higgs large 6± is needed, while 6a could 
be small. In this work we limit ourselves to the light SM-like Higgs boson hs- 

LEP II limits As D^, Da, Dh do not couple to fermions, the LEP limits based on Yukawa in- 
teraction for the standard 2HDM don't apply. However, the signatures are similar to neutralinos 
and charginos interactions in MSSM. The absence of a signal within searches for supersymmetric 
neutralinos at LEP II was interpreted within the IDM in paper [T7]. This analysis excludes the 
following region of masses: Md^ < 80 GeV, M^^ < 100 GeV and ^a > 8 GeV. For 5^ < 8 GeV 
the LEP I limit Md„ + Mdj, > Mz applies. 

3 DM relic density constraints 

In this analysis we assume that Dh is a dominant component of the observed DM, with its 
energy relic density in the Universe estimated to |19) : 

riDMh^ = 0.112 ±0.009. (13) 

Various studies show [TT| [T21 [131 [HI [13 [IS] that for IDM in most regions of the parameter 
space VtoMh'^ is too low to fulfil the astrophysical constraints. However, there are three allowed 
regions of Md^: (i) light DM particles with mass close to and below 10 GeV, (ii) medium DM 
mass of 40 — 80 GeV and (iii) heavy DM of mass larger than 500 GeV. 

Those regions are further constrained by the value of A345, which governs the main decay 
channels (decay through Higgs exchange) for M^^ < My/. In general, for larger IA345I the relic 
density decreases due to the enhanced D^Dh annihilation via hs into pair of fermions (typically 
hh). 

On the contrary, the value of A2 does not influence the DM relic density directly. At the 
same time, this parameter is difficult to access at colliders. Therefore, this parameter is usually 
fixed to abitrary small value in the DM analysis of IDM [TTl \n\ [T3| IHl [TSl [IS]- However, as we 
argue in this paper, value of A2 limits the value of A345 by the positivity constraints. Its value 
is also limited by the constraints arising from the possible existence of extrema and vacua of 
different properties than Ji during the evolution of the Universe. Therefore, in this analysis we 
will consider the constraints for Mojj as well as for both A345 and A2 arising from the astophysical 
data. 

In this work we derive the corresponding relic-density constraints for IDM, using micrOMEGAs 
|20) . with the IDM implemented by us. In analysis we respect all other existing limits and confirm 
findings of [TTl [T2l [T3l [HI [IHl [IS]. The micrOMEGAs program neglects temperature dependence 
of physical parameters and a possibility of more than one phase transition. In paper [6j we 
concluded that if in the past there were sequences of phase transitions, then the Universe en- 
tered the inert phase with DM candidate at lower temperatures that in the typical one-stage 
EWSB. In principle, this should be considered while solving the Boltzmann equations for DM 
relic density, especially while the final phase transition into /i takes place at the low tempera- 
tures during the freeze-out. Also the latent heat of the Ist-order transition may give significant 



corrections in some regions of the allowed parameter space [21j . In this sense, the energy relic 
density calculations in this paper should be considered only as a preliminary estimate]^ 

Below we discuss the general trends in the change of Vlj^Mh"^ value in three regions of DM 
masses. In all examples we consider the light SM-like Higgs with its mass fixed to Mh = 120 GeV. 
We expect the VLohih"^ to be in the 3cr WMAP range, namely: 

0.085 < VtDMh'^ < 0.139. (14) 

Note, that the medium DM mass region is strongly constrained by the existing collider data. 



3.1 Low DM mass region 

In the low mass region M^^ is of the order (4 — 8) GeV and masses of D^ and D^ are almost 
degenerate with 5a ~ 5h± ~ 100 GeV [llj . 

Large mass splittings between the Dh and other scalar particles do not allow for the coan- 
nihilation. The main decay channel is D^Dh -^ ff {cc pair for M/j^ = 4 GeV and hh pair for 
higher M^^) via Higgs boson exchange. Small IA345I generally gives high Qj^Mh'^ well above the 
WMAP limit, as this decay channel is suppressed. The larger IA345I is, the lower Q^Mh"^ gets. 
The WMAP allowed region of IA345I is around (0.4, 1.2), but this strongly depends on the exact 
value of Mdjj (figure [I|. Note, that the 1 GeV difference in mass is this region causes significant 
change in the VLouh'^ for a given A345. 

One should remember that changing of A345 is allowed only in the region allowed by the value 
of A2 (in the calculation of Q^iMh"^ we set A2 = 5 which corresponds to 1.5 > A345 > — 1.5). In 
general, as we allow the higher values of A2 the range of A345, over which we can scan, extends 
and for a chosen M^h the higher IA345I is the lower QoMh'^ we get. 

For the low mass region relic density does not depend on the value of masses of the much 
heavier Da and D^. In our analysis we keep the mass splittings fixed to 6a = 100, 6± = 105 GeV 
and M^^ changes in the allowed region (4 — 8) GeV. As M^^ decreases, the value of relic density 
for the chosen A345 grows, as shown in the figure [T} 
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Figure 1: Relic density for M^^ = (4-8) GeV, A2 
Horizontal lines denote WMAP 3a allowed region. 
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100 GeV, ^4 



105 GeV. 



^More work is in progress. 



3.2 Medium DM mass region 



In the medium mass region the mass of Dn is in range M^j^ = (30 — 80) GeV. Mass sphttings 
5 A-, 5± can be of the same order 5^ ~ (5± = (50 — 90) GeV, but also small values of 5a are possible 
(of the order of 10 GeV) ^. 

The medium mass region Qoivih"^ is very sensitive to the value of Md^. Therefore in this case 



it is difficult to make a general statement similar to the low mass case (section 3.1). However, 
we see some regularities in this behaviour, related mostly with the effects of coannihilation, as 
discussed below. For our calculation we use A2 = 0.3, which in this mass range allows scanning 
over A345 in range ~ (—0.3, —0.35). 



Small 5a Let us first consider the small mass splitting between Dh and Da (figure 2a): 

5a = 8 GeV , 5h± = 50 GeV . 

Small 5a makes the coannihilation {Dh, Da) important, decreasing the Vtohih'^ below the WMAP 
limit in most of the A345 parameter space. Usually only the region close to A345 = is inside the 
WAMP limit. For the larger masses of Dh-, larger A345 are allowed. 

Note, that the smaller masses of D^ which may give a different results are heavily constrained 
by the LEP 11 data. 

Large 5a If mass splittings between Dh and other scalars are larger {5 ~ 50 GeV) then the 
coannihilation is no longer important. For smaller masses and large mass splittings the behaviour 
is similar to the one in the low mass region (A345 region around is excluded and larger IA345I 



are allowed, figure 2b). The allowed values of IA345I are larger than in the case of small 5a- 



This region of masses is the most natural to consider as the constrains for the masses of 
scalars are not so tight as in other cases. 
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Figure 2: Relic density for medium mass region, results for A2 = 0.3. 
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3.3 High DM mass region 

In the high mass region Md^ = (500 — 1000) GeV all dark scalars have almost degenerate masses 
and the mass splittings are small of the order oi 5 ^ 5a ^ 5± < 12 GeV due to the perturbativity 
conditions [TT] . 

In this high mass region, with such small mass splitting, the coannhilation between all dark 
scalars particles is very important. The results show very high sensitivity to the value of (5's, as 
shown in figure [3} Here we fix M^h = 800 GeV, while 5 varies from 1 to 10 GeV, respectively. 
In general, as 5 grows, the value of relic density for given A345 decreases. For 5 = 10 GeV, for 
every value of A345, we are below the WMAP limit. 
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Figure 3: Relic density for high mass region: (A345, QoMh"^) plot for fixed DM mass M^)^ 
GeV and 6a^± = (1, 3, 5, 7, 10) GeV, resuhs for A2 = 1.5. 
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In figure ffl Md^^ changes between 500 and 1000 GeV, while 5's are fixed. We first consider 
6 = 1 GeV (figure 4a), where lower values of Md^ correspond to the lower values of Q^Mh'^- If 
5 = 10 GeV (figure 4b) then for most of the chosen values of Md^ we are below the WMAP 
limit. Mass around 1000 GeV gives the proper relic density for higher values of IA345I, note that 
A345 < 0. 
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Figure 4: Relic density for high mass region for varied DM mass. Results for A2 = 1.5. 



4 Thermal evolution 



4.1 Thermal evolution of the Universe 

Evolution of the potential Following the approach presented in O HI [5l |6l [7] we consider 
the first nontrivial temperature corrections to the potential V. The coefficients X's of the quartic 
terms are unchanged, while the quadratic terms vary with temperature T with coefficients: 



ml{T) 



m 



11 



CiT\ mUT) 



C\ 



3Ai 



2A3 + A4 ^ 3^^ + g 



/2 



9l^gl 



C2 



"^22 - C2T^ , 
3A2 + 2A3 + A4 



+ 



3^2 + g' 



(15) 



6 8 2 ' ^ 6 

The EW gauge couplings g and g' and Yukawa couplings for h and t quarks are: 



g = 2Mw/v = 0.652, g' = 0.351; gi = V2mi/v {gt ^ 0.99, gt ^ 0.02). (16) 

As shown in [6j, depending on the value of X's each of the coefficients Ci and C2 can be either 
positive or negative. In this work we don't discuss the possibility of non-restoration of the EW 
symmetry in the past [22] and therefore we consider C2, Ci > only |6]. 



Extrema and vacua 


name of extremum 


vev's 


properties of vacuum 


EW symmetric: EWs 


vd = 0, vs = 


Massless fermions and bosons and massive 
scalar doublets. 


inert: h 


n 2 2 ^^11 


Massive fermions and gauge bosons; scalar 
sector: SM-like Higgs hg and dark scalars 
Dh,Da, D^ with DM candidate Dh- 


inertlike: I2 


n 2 2 ^^22 
^2 


Massless fermions and massive gauge 
bosons; scalar sector: Higgs particle 
h£, (no interaction with fermions), four 
scalars Sh-, Sa-, S^ , no DM candidate. 


mixed: M 


vl,vl>0, v^=vl + vl 


Massive fermions and bosons, 5 Higgs par- 
ticles: CP-even h and H, CP-odd A and 
charged H^ , no DM candidate. 



Table 1: Extrema and properties of the vacua [6J. 

Extrema during the evolution Depending on the values of the parameters the D-symmetric 
potential (tsl) may have different types of the vacuum state. The general form of an neutral 
extremum is given by the following solution of the extremum conditions: 



{^s) 
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{^d) 



y/2\VD 



{v^ 



+ vh). 



(17) 



The properties of these extrema, provided they are realized as vacua, are listed in the table 
[l| As we discussed in ^, each of this vacua can be realized in the separate regions in (/ii,/i2) 
plane. 

Types of evolution As the quadratic terms in V vary with T, the vacuum state changes. As 
time grows (and temperature decreases) the vacuum state of the Universe changes along one of 
rays in the phase diagrams shown in figure [5} where the allowed regions for all vacua EWs, Ji, 
I2 and M are shown as well. 

The types of evolution (rays), along with conditions on R and fii, ^2 that need to be satisfied 
in order to have a chosen type of evolution, are summarised in the table [2] in the appendix. 
Note, that in existing in literature DM analysis only the rays that correspond to a single phase 
transition EWs — )■ Ji have been considered. 



4.2 Rays in the (A345, A2) plane 

To get a deeper insight into evolution of the Universe we consider now parametrization of rays 
by the physical parameters like masses and self-couplings. We will use the (A345,A2) plane, to 
represent regions where a certain ray is realized, treating masses of scalars M^^, M^^, M^^, M£)± 
as the input parameters. Conditions from table [2] have been translated into conditions for 
(A345,A2) parameters and are presented in table Islin the appendix. Note, that only for the 
case of a single phase transition we need to impose the constraints for A345 parameter. In case 



of sequences with two or three phase transitions, constraints on the quartic couphng A2 are 
sufficient. 

The (A345, A2) phase space can be divided into separate sectors. Let us introduce A345 param- 
eter, which separates these sectors: 

A,,., = —^. 18 



V 



For A345 < A345, which corresponds to ^2 < 0, the only EWv extremum, and thus the vacuum, 
that ever existed is the inert phase Ii. The only possible sequence is EWs -^ Ii. For A345 > A345 
(/i2 > 0) other EWv extrema (M, I2) could be realized in the past. Their realization as vacua 
depends on the value of A2. 

Distinct regions of A2 are separated by A2, which depends on A345 and other parameters as 
follows: 

~ _ -jSg^ + g") (MB, - 3M|,^ + 2M^ J ^ A345 , (2M|,^ - t;%45) 
^' - 4 3^^ +^+ 12 "" 

{-4Ml^ + 12MI^ - ml^ - {9g^ + ^'^ + I2gl + 12^,^ + 8A345)t;^) 

v^l ■ ^'^^ 

Note, that for A2 > A2 only one phase transition may be realized, while for A2 < A2 sequences 
with two or three phase transitions are possible. 

Below we perform a general analysis of possible rays on the (A345, A2) plane for various masses, 
assuming that today Ji vacuum is realized. 

4.3 Benchmark points 

For further discussions we will use the reference points Bl — BA chosen in three regions of Mjj^ 
as follows. For all points Mhg = 120 GeV. 

• Low DM mass region: for the reference point Bl we have: 

Bl : Md^ = 8 GeV, 6a = 100 GeV, 6± = 105 GeV. (20) 

• Medium DM mass region: as discussed, large or small Sa are possible. Those are the 
benchmark points B2 and B3, respectively. 

B2 : Md^ = 45 GeV, 6a = 70 GeV, 6h± = 70 GeV, (21) 

B3 : Mdh = 70 GeV, 6a = 8 GeV, 6h± = 50 GeV. (22) 

• The high DM mass region set BA is defined as: 

BA : Mdh = 800 GeV, 6a = 1 GeV, 6h± = 1 GeV. (23) 



Figures (pa-d) show the regions, where concrete rays can be realized in (A345, A2) plane for 
mass points Bl —BAq In general, the higher Djj mass is, the higher values of A345 are needed for 
the realization of the sequences with more than one phase transition. More discussion is given 
below (sec. 4.4 4.6[), where also the relic density as the function of A345 is discussed. 



^Note, that the scales in those figures are different. 



4.4 Low DM mass region 



Rays in (A345, A2) plane Figure 6a sliows tlie (A345, A2) plane for the low mass point Bl. In 
this case possible types of evolution are very limited. There is an allowed region for the single 
phase transition (realized by ray Ic for R < and ray lib for 1 > i? > 0) and sequence of three 
phase transitions with transition through M (ray VI). Important feature of this low DM mass 
region is the fact that Ist-order phase transition EWs -^ I2 -^ h is not possible (rays IV and 
V), as the allowed phase space for those rays shrinks to the degenerate line A2 ~ X^^^v'^/Ml . In 
short, fulfilling R> 1 condition is not possible for the low DM mass region. 



Relic density The obtained energy relic density as a function of A345 is presented in figure 7a 
for a low mass benchmark point Bl. There are two regions in agreement with the WMAP limit: 
A345 e (—0.405, —0.315) and A345 e (0.315, 0.405), presented as the dark regions in (A2, A345) plot 
(figure [7b) ). For the negative values of A345 ray Ic, for positive - ray VI or ray lib, depending 
on the value of A2, can be realized. Note, that the value of A2 should be quite large to have 
any kind of solution that fits into the WMAP region. Because of the positivity constraints 
for A345 G (—0.405, —0.315) we need A2 > 0.45 if we want to have ray Icrj For A345 > and 
1.9 > A2 > 0.45 only ray VI is realized; if we want to have ray lib we need A2 > 1.9. 

Calculation of Q^Mh"^ was done for A2 = 5, which does not influence its value. However, 
as we discussed above, the value of A2 is limited by the other constrains for certain rays. It is 
limited by the positivity constraints A345 + v^AiA2 > for A345 < 0, while for A345 > positivity 
constraints don't give the further constraints on A2, but its value is limited by the conditions of 
realization of particular rays (table |3]). 

4.5 Medium mass region 



Rays in (A345,A2) plane Figure 6b shows the allowed rays for the reference point B2 in the 
medium mass region for Md^ = 45 GeV. In this case all types of evolution are possible. Note, 
that rays that correspond to the sequence of two or more phase transitions (rays IV, V and VI) 
appear for larger values of A345. For the point B3, with Md^ = 70 GeV, the types of rays are 



similar to the B2 case, as shown in figure 6c Some small differences are visible, for example 
ray VI requires larger value of A2. Those differences are mainly due to the value of M/)^, mass 
splittings 6a, S± are less important (see section 5). 

Relic density Point B2 is similar to the low mass region as here the coannihilation is not 
relevant. However, due to the larger Dh mass the annihilation through Higgs exchange is more 
efficient and the obtained A345 values are smaller. Again there are two allowed regions that give 
the proper relic density (figure [s]). Note, that this time for the positive A345 rays lib. Ha and HI 
(the later with the coexisting local minimum I2) are allowed by the relic density data. 

In case of the point B3 due to the small value of 6a the coannihilation {Dfj, Da) is important 
for the relic density. Its effects are clearly visible for the very small values of A345, where 



annihilation through Higgs is strongly suppressed (figure 9a). We get the proper relic density 



for A345 G (—0.1,0.1) with the small region around excluded. This corresponds to rays la-c. 



depending on the value of A2 (figure 9b). Here, sequences with more than one phase transition 
give ilDAjh^ below the WMAP limit. 



^For positive allowed A345 and A2 < 0.45 the Universe is in the state of intertlike vacuum. 
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4.6 High mass region 

Rays in (A345, A2) plane For small values of A345 only rays Ic, lb and la are possible (figure 



6d). The types of evolution other than a single phase transition could happen only for the very 



large values of A345 > ^Ml^/v"^, which for Md^ = 0(800) GeV gives A345 = 0(20). 

Relic density In high DM mass region, due to the very small mass splittings coannihilation 
processes are the most important ones. For chosen set of masses the allowed ranges of A345 are 



0.42, —0.28), (0.21,0.35), as shown in figure 10a It is clear that for this set the allowed rays 



are la, lb, Ic (single phase transitions). The other rays could be realized for A3 > 22, but first 
of all, it would violate the perturbativity constraints. On the other hand, the value of Q^Mh"^ is 
very low (of the order of 10^'^) for such large A3 (and A245). 

5 Sensivity of the types of evolution to the dark scalar 
masses 

The realization of different types of evolution depends strongly on values of all parameters in 
the potential, as presented in numerical examples in |7j. However, in general the value of Md^ 
affects the type of evolution stronger than the values of M^^ and Mij± . 

Let us consider first the values of masses for dark scalars allowed by the relic density con- 
straints for the low and medium mass of M^^ ^ (8 — 80) GeV: 

Md^ = 110 GeV, Md± = 110 GeV. (24) 

Figures (|TTk-d) show the regions of parameters on (A345, A2) for the different values of M^^ 
(from 8 GeV to 60 GeV). For the low mass region single phase transition (rays Ic, lib) and 
transition through M vacuum (ray VI) are possible, other types of evolution appear for the 
higher values of M^,^. As the mass grows, the region of ray 11a expands and thus ray VI 
becomes possible for the larger values of A2. 



If we fix the value of Md^ and let the other masses vary (figures 13i-d) the picture is not 



significantly affected. This is presented in the example of the large mass splittings Sa, 5±, however 



using small 6a doesn't make much difference (compare with figure 14a noticing the difference in 
scales) . 

Low DM mass region Since M^h is small compared to the other masses, i? > 1 is in fact 
not possible and the possibilities presented in table [s] are reduced to only three rays (Ic, 11a, 
VI). In the considered range of M^u the main mass parameter that influences the boundaries of 
rays is mass of the Higgs, Mh^. Changing M^^ in the allowed range with fixed Sa,S± for this 



region does not change the picture [6aj Boundary line between rays lib and VI depends mostly 
on Mhg and fixed gauge and Yukawa couplings. Small corrections (AA345 = 0.01) coming from 
the change of 6a and 6± arise for the small (A2 € (0.01, 0.1)) and large (A2 > 5) values of A2 . 

In section 3 we presented the fi^ijv/^^ calculations for this region. Here we see that the values 
of Mjj^ and A2 are the most important parameters. In a way, A2 is more important than A345, 
because even if we are in the allowed region of A345, the small value of A2 may lead us to a region 
excluded by the positivity constraints. As Md^ grows, we enter the WMAP region for the lower 



absolute values of A345 which corresponds to lower values of A2 (figure 12). 
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Medium DM mass region For large mass splittings 5a ~ 5± all rays (and all types of 
evolution) are possible in the range M/)^ = (30 — 60) GeV (figures 6b|llb-d). As mass grows, 



ray VI requires higher values of A2 and the region of this ray shrinks in benefit of rays Ila, IV 
and V. For the medium mass region the change of Mojj effects the evolution for low values of 



A2 more than the change of M/j^ and M^ii for fixed value of Mr,^ (figures 13i-d). 

For small 5a the picture is similar to the previous one, however as the M/?^ grows the allowed 
region of ray VI becomes limited in benefit of rays Ila and IV. For M^i^ = 80 GeV ray VI is still 



possible but it requires higher values of A2 ~ 0.7 (figures 14a,-c). 

The types of rays that may be realized with the A345 inside the WMAP limit depend heavily 
on the specific value of Md^- 

For small 5a region around A345 = gives proper VtoMh"^- Here only sequences with a single 
phase transition can be realized. As M/j^ grows, larger A345 are allowed by WMAP. However, 
also the region where the more complex sequences can be realized shifts towards even larger A345 



(figure 14). In general, this case may provide us only the EWs — )■ Ji sequences in agreement 



with WMAP (figure 15). 



For large 5a also sequences HI- VI can be realized with the proper relic density (figure 16). 
This is the case for masses Ma^ ^ 40 GeV, where rays with the Ist-order phase transition are 
possible. For this value of mass the WMAP allowed range is in A345 G (0.13,0.17). One should 
keep in mind that this region can be constrained soon by the DM direct detection data. For 
larger masses the region of those complex sequences is shifted towards the larger values of A345. 
In general they give VL£,Mh? below the WMAP limit due to the DhDh -^ hs -^hh channel. In 
this case we are left only with the EWs — ?■ Ji sequences realized for smaller A345 inside WMAP 
limit. 

High DM mass region In this region the lower value of M/)^ makes the realization of other 
than a single-phase transition types evolution possible for lower values of A345. Still, for M/j^ = 
500 GeV this value is high of the order 0(10). If we limit ourselves to Aj ~ 0(1) then only single 
phase transition is possible in the high mass region and the (A345, A2) plot is virtually identical to 



figure 6d As discussed before, proper VtoMh is for A345 < 0(1) and thus the complex sequences 



are not in agreement with WMAP. 

6 Conclusions and outlook 

In this work we analysed the possible types of evolution of the Universe (rays) in the IDM, 
including constraints from the DM relic density data (WMAP). Three regions of the allowed 
dark matter masses: low, medium and high mass were considered. Evolution wasparametrized 
by masses of scalars (including dark matter particle) and two parameters (self-couplings) A345 
and A2. With fixed scalar masses we can find the values of (A345, A2) that satisfy the conditions 
of a chosen ray (table [3] in the appendix). We have shown that those three regions of DM mass 
exhibit the different behaviour: both in the possible types of evolution and the energy relic 
density values. 

For the low DM mass region only three kinds of rays can be realized: ray Ic, lib and VI. 
This is due to the very low Dh mass and the exact values of the much heavier other scalars do 
not affect the picture. The energy relic density also doesn't depend on M^^ or M£i±, but it is 
stronly influenced by the exact value of Md^. In general, rather larger IA345I ~ (0.4 — 1) are 
needed to fulfill the WMAP constraints. Due to this fact, in our approximation of the Q^Mh"^ 
calculation, sequence EWs — )■ /2 — ?■ M — )■ Ji gives the proper QuMh"^- 
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For the medium DM mass region all types of evolutions are possible. Again the main param- 
eter that affects the regions in (A345, A2) plane is the mass of DM particle, however for lower A2 
also the effect of changing of M^,^, Md± are visible. In this region Q^Mh'^ depends heavily on 
the value of Md^, but also on Sa- There are two separate regions. For small 6a coannihilation 
is possible and in general small IA345I give the proper DM relic density. Large 6a case is similar 
to the low DM mass case as the coannihilation is not important and the allowed IA345I tends to 
shift towards the larger values. 

In the high DM mass region all dark scalars have almost degenerate and very large masses. 
This has two main consequences. First, the (A345, A2) plane for A < 0(1) is reduced only to the 
regions corresponding to the simple sequences with a single phase transition. Second, the energy 
relic density is strongly lowered for most values of IA345I, because of the coannihilation. 

We would like to comment, that the high mass region requires very careful fixing of the 
parameters, as it is very sensitive to changing the values of 5's. This certainly doesn't seem to 
be natural (fine-tuning) . From this point of view the most interesting is the medium mass region 
with the large mass splittings, as here the adjustment of the parameters isn't so strict. 

It is worth to stress the importance of the often neglected A2 self-coupling. We argue, that the 
astrophysical data should be used to limit the values of this self-coupling. Firstly, the Q^Mh"^ 
calculation for a fixed M/j^ gives the allowed values of the A345 parameter. This calculation 
doesn't depend on the exact value of A2. Secondly, those obained values of A345 should be 
used to constrain the A2 parameter through the positivity constraints or the conditions for the 
realization of the different rays. Therefore, fixing A2 to an arbitrary value during the calculations 
may result in the exclusion of the WMAP allowed region. This is especially visible in the low and 
high DM mass region, where for the A345 < the positivity constraints require A2 > (0.2 — 0.6). 

As we mentioned in section 3, the calculations of Vt^Mh'^ in this work should be treated 
only as an estimate. We do not take into account the corrections to those calculations arising 
from the thermal evolution of the Universe, such as the infiuence of complex sequences of phase 
transitions or the effect of the 1-st order phase transitions. Also, in certain ranges of (A345, A2) 
for low and medium DM mass region the final phase transition into inert vacuum may take 
place at the lower temperatures, which may affect the freeze-out. Those effects are now under 
investigation. 

Acknowledgments I would like to thank I. Gizburg, K. Kanishev and M. Krawczyk for co- 
operation and discussions. I would like to thank M. Krawczyk for reading the manuscript. Work 
was partly supported by Polish Ministry of Science and Higher Education Grant N N202 230337. 

A Conditions for rays 

Below we present the tables that contain the conditions for the realization of different types of 
evolution presented in two planes: [fii, H2) in table 2 and (A345,A2) in table 3. Here we use the 
redefined evolution coefficient: 



c = C2/C1 = (c2VAi)/(ciVA2). (25) 

To make the table more transparent we use the previously introduced boundary values A2, A345: 

_ -jSg' + g'') jMl^ - 3Ml^ + 2Ml^) ^ A345 , (2M|,^ - t;%45) 
^' - 4 3^^ +^+ 12 "" 
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Ray no. 


R 


Conditions 


EWs -^ h 


la 


R> 1 


/i2 <0 


Ila 


R> 1 


< /i2 < Mm{fiiC,fiiR-^) 


III 


R> 1 


fiiR~^ < yU2 < Min(/iic,/ii) 


lb 


1 >i?>0 


/i2 <0 


lib 


1 >i?>0 


< /i2 < Mm {fiic,fiiR) 


Ic 


0>i?> -1 


fi2 < fiiR < 


EWs ^h^ h 


IV 


R> 1 


yUlC < yU2 < yUl^ ^ 


V 


R> 1 


Max (/iic, jiiR^^) < /U2 < /ii 


EWs ^h^ M -^ h 


VI 


1 >i?>0 


/ilC < /i2 < yUli? 



and the abbreviation: 



Table 2: Possible rays ; (/ii,/i2) conditions. 



A 



2M 



Df 



345 



a 



(t^^A345 



2m: 



i^f 



Ml 



(27) 
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Figure 5: Possible rays on (yUi,/i2) plane. Vacua: EW symmetric (white region), inert (yellow, 
lightest shade), inertlike (blue, dark shade) and mixed (green, medium shade). Vertical hatch 
between lines A and B in figure |5a] denotes the area when I\ is a local minimum, while I2 is the 
global one; horizontal hatch between B and C - /2 is a local minimum, while I\ is the vacuum. 
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Figure 9: Relic density for B3 with Md^ = 70 GeV, 6a = 8 GeV, 6± = 50 GeV. 
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Figure 12: (A345, A2) plot for low DM mass region with the WMAP-allowed areas for Md^ 
(8,5,4) GeV, respectively. 
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